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Abstract 

Effect of disorder in metallic thin film is examined as a possible mechanism of 
the Superconductor-Insulator (S-I) transition. The critical value of disorder corre- 
sponding to the transition point is found analytically by using Matsubara-Matsuda 
model and Green's function method. 

1. INTRODUCTION 

At zero temperature two-dimensional systems of interacting electrons are speculated to 
show a quantum phase transition between superconducting and insulating (S-I) phases. 
S-I transition can be driven by tuning some parameters such as the disorderS^i, carrier 
concentration!'!, magnetic fieldi'l, dissipation!!!, and so on. It has been observed in granular 
and amorphous films such as Bi, Pb, Josephson junction arraysIH0, and 4 He absorbed 

on porous media0, in two dimensions. Here especially we are interested in superconducting 
films. 

Superconducting films are made by repeated small increments of materials onto substrates 
held at low temperatures in an ultra-high vacuum. Then the system is found to form super- 
conducting islands separated by thin insulating regions. The films are strongly disordered 
due to distribution of the island size and the coupling between islands. We assume, in this 
paper, that the randomness in island potential is crucial in determination whether the film 



is superconducting or insulating. 

Natural questions arises : how much disorder is necessary to destroy the superconducting 
properties? Precise answer to this question will depend on details of the interaction among 
the particles, as well as the nature of the disorder in the system. 

Destruction of superconductivity with disorder can be caused by localization effect rather 
than Cooper pair braking. Then the system becomes insulator on the macroscopic scale in 
spite of local superconducting correlations. Relying upon a universality hypothesis, we 
regard Cooper pair as a tightly bound hard-core boson as long as we restrict ourselves to 
the vicinity of S-I transition point, disregarding any microscopic detail such as electrons, 
phonons, and their interactions. After all our problem becomes equivalent to dirty boson 
problem!0ll, which has been extensively studied using quantum Monte Carlo simulationsSi, 
real-space renormalization-group techniquesBIll, strong-coupling expansions^, and other 
ways. 

In this paper we shall adopt Matsubara-Matsuda model0"0 which was established for 
tightly bound hard-core bosons, and calculate the critical value of disorder which corresponds 
to the transition point with double time Green's function&~il. 



2. MODEL 

Inspired by the work of Matsubara and Matsuda (M-M) who formulated the model of 
hard-core boson in lattice era we introduce Bose-Hubbard type Hamiltonian : 

H = U^2 ^ifij - t^2(a\a i+1 + diaj + i) - XX ^ + $h)hh (!) 

i,j i i 

where aj, fij, and hi are the creation, annihilation, and number operators of Cooper pairs at 
the z-th lattice point, respectively. We exclude boson double occupancy by assuming infinite 
repulsive on-site interaction. We have introduced a finite repulsive interaction U, for two 
nearest-neighbor bosons. In addition, we restrict hopping t to nearest-neighbor only. In this 
model we introduce random site potential Shi which represents disorder due to distribution 
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of island size. It takes both positive and negative value. The external potential is described 
as the sum of uniform part h and random part Shi. The boson occupancy is restricted to 
one or zero. In other words each site has only two states, empty or full. The hardcore boson 
system is therefore isomorphic to a spin 1/2 system. Each operator of Cooper pairs can be 
replaced as spin operator!! : 

aj = Sr = i(fif - zaf), @a, = 5+ = |(af + to?), (2) 

fH = \-St = \{l-atl (3) 

where af, of, and of are Pauli matrices for the 1/2 spin at the z-th lattice point, respectively. 
From (|l]), (@), and (|3|), the Hamiltonian becomes 

J z = -4C7, @ = 2t, (5) 

hi = h + Shi. (6) 

As a result, interaction energy and kinetic energy are now represented by spin exchange en- 
ergy, and chemical potential is represented by Zeeman energy. Our system can be described 
with an anisotropic ferromagnetic Heisenberg Hamiltonian. 

Therefore we can interpret ordering of pseudo-spins in x-y plane as superconductivity and 
z direction as local density nuctuationlll. 

(at), (a) & (S-),(S + ), (7) 

(n) ^ (S z ). (8) 

Since (S + ) and (S~) represent superconducting order, their value can be regarded as a 
criterion to judge whether the system is, superconducting or insulating. So we shall calculate 
them in the next section. 



3. CALCULATION 



The problem now involves statistical physics of the pseudo-spin system described by the 
ferromagnetic Heisenberg model. The temperature dependent retarded Green 's functionEM 
with two operators S^, Sr and a coefficient hi, is introduced as 

«^(t);^T(0)» = -ieWdhiStMJriO)]). (9) 

(t) is a Heisenberg operator at time t ; 6{t) is the step function ; square brackets [■ ■ •] 
denote commutators, and single angular brackets < ■ ■ ■ > denote thermal averages. A 
straightforward calculation yields the time-Fourier transformed equation of motion for the 
Green's function (9) 

cudhS^Sr)) = (M+ £7]) + <<M+ W];£r>>. (10) 

In our two-dimensional square lattice model, we only consider coupling between nearest- 
neighbor sites (see Fig.l), i.e. site % with site % + a n (n = 1 ~ 4). 

Sr)) = ([htSf, Sr}) + S-)) 

+ J2(U(hiStS* +a ;S.)) -J^dhStSt^Sj))). (11) 

a 

We can solve the equation of motion ( |TTD easily by recovering translation invariance 
by impurity averaging procedure, and this can be done by averaging the magnetic field as 
follows 

{h i ) = (h + 8h i ) = h, (12) 
(hf) = (h 2 + 2h5h l + 5h 2 ) = h 2 + 5h 2 . (13) 

Since 5hi takes both positive and negative values, the first order average equals to zero. 
Then the equation with translation invariance and randomness contribution is obtained as 

cuh((S+; Sr)) = h([S+, Sr]) + (h 2 + 5h 2 )((S+; Sr)) 

+ £ (j z h((StS? + a, ^)) - J xy h((S-Sf +a , Sr))) . (14) 



For simplicity, we decouple higher order Green's function in the following manner 



({Si S* +a ; Sj)) — > (S z )((St; S*)), 



(15) 



<<^+ a ;^>) - (S z )((S+ a ,S*)). 



(16) 



Then (|14) becomes 



w«fl?-; 1 ffir» = 2<5*> + 



/i 2 + Sh 2 
h 



((St; Si)) 



+ E (US z )((S+;Sr)) - ^(^)((^ +a ;^-))) . 

a 

Fourier transformation ( |i~7[ ) with respect to space results 

"((st; sz k )) = 2(s*) + U—((s+ ] sz k )) 



h 



Y.^ ka )((st;sz k )). 

a J 



(17) 



+4(S Z ) \j z -^L^r~**\ //<?+■ 

Putting them in order, we obtain the formula for the Green's function ((S£; SZ k )) 

2(S Z ) 



(18) 



G(k,u) = ((S+;SZ k )) 



where 



lu - A(S z )uo k - H 



(19) 



Aka 



(20) 



H 



h 2 + 5h 2 
h ' 



(21) 



Since the poles of Green's function indicates the excitations of the ferromagnetic spin wave, 
we see that two-body interaction J z , external field h, and Ah assist crystallization and 
hopping integral J xy causes quantum fluctuation. 

Next we treat the correlation function (Sj Sf). It can be calculated by using the spectral 
theorem 



(S~S + ) 



1 f w 
4^ U dK 



dk v — / duo 



{<%{<*)- Gt{u,)} 



er — 1 



(22) 
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where Gf(u), G k (uj) are retarded and advanced Green's functions. From equation (|T^) we 
obtain 



G%(w) - G£(w) = -2m6(u - A(S z )u k - H), (23) 
then the integration of the right-hand side of with respect to uo becomes 

- - (S Z ) = -P-r dk x dk y * . (24) 

2 47r 2 i-. X J-* y exp[{A(S z )u k + H)/{2T)]-l 

Putting them in order, we obtain the relation 

^ r kdkcoth + (25) 



71 Jo V 2T 

Since low temperature region is of our main interest, we neglect large spin wave vector. This 
approximation gives uo k as 

u k - ^Jxyk 2 + J z - J xy . (26) 
Therefore the excitation energy gives 

u = (S z )J xy k 2 + 4(S Z )(J Z - J xy ) + H. (27) 

The second and third term on the right-hand side means the gap of the spin wave excitation 
A at k ~ 0, which is caused by anisotropy and external field. Here two conditions must 
be fullfilled in order that the ground state is superconducting. The first one is that A is 
sufficiently small compared to T, i.e. 

A = 4(S Z )(J Z - J xy ) + H ^T, (28) 
the second one is 

Jxy ^ Jz- 

(29) 

Using the above conditions, the integration on the right-hand side of (^Sj) can be calculated 
analytically, and so we obtain the expression at low temperature 
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(s z ) 



i7( e -Jx y /T _ X ) 



H 



ir 2 J xy -4(J z -J xy )(e^y/ T -l) 



4(J, 



<7rj/) 



(30) 



Due to rotation symetry of spin in x-y plane, we can select the X~cLXlS clS the grand state di- 
rection without lost of generality. Then we obtain the relation between the superconducting 
order parameter (S x ) and randomness Ah 



(S x ) 



r 2 

2, 



(fry 



i 



H 2 



2 ^ ' 1 ( '^2 Jxy) 



(31) 



4. RESULT 



Before examining the physical meanings of fl3~T|), we introduce some normalized factor 



rj = 



Ah 



m = 



h 



(32) 



'xi/ u xy <Jxy 

Which corresponds to normalized strength of randomness, two body interaction, and chem- 
ical potential, respectively. Then ([H]) becomes 



(S x ) 



1 - 



(e + v 2 ) 2 



(33) 



2^ 4£ 2 (l-/t) 2 ' 

This equation describes the behavior of the superconducting order parameter (S x ) with 
respect to the strength of randomness 77 (see Fig. 2). From Fig. 2 we can see that (S x ) is 
damped by increasing rj. By expanding (|3"3]) in power of rj 



2\ 4(1 -k) 



77 



4(1 -k) 2 -f 2 



(34) 



we can see the damping of (S x ) occurs proportional to rf. At a specific point, (S x ) completely 
disappears, which implies the existence of a phase transition. As mentioned in sec .2, we 
interpret phases (S x ) 7^ as superconducting and (S x ) = as insulating. Here the central 
result is the existence of the critical value rj c for the strength of randomness. r] c divides the 
two phases (S and I) and gives the transition point. From (|33"D rj c is given by 



2£(i 



(35) 



From Fig. 2 we can see that the strength of two-body interaction k advances the transition. 
Therefore k can be also regarded as a parameter of S-I transition, (see also Fig.3). 

Varying k instead of 77 (see Fig.3) also gives a phase transition, this time it corresponds to 
a Mott transition. Especially when 77 = 0, (S x ) sharply decreases at k ~ 1. This corresponds 
to a transition in an 2D isotropic system i.e. 77 = 0, k = 1. This is based on the general fact 
that (S x ) doesn't exist in two-dimensional system. 

Now we vary both r\ and k as a parameter of the phase transition, the border line between 
two phases (S and I) is determined by the condition that superconductivity disappears (i.e. 
(5*) = 0). 

From (|33| ) this relation corresponds to 

1 - .%, V \ 2 = 0. (36) 
4£ 2 (1 - k) 2 

Real (S x ) corresponds to the system being superconducting, in this case 77, k have to satisfy 



(£ 2 + ri 2 ) 2 



On the other hand, non-real (S x ) corresponds to the system being insulating, this time rj, k 
have to satisfy 

(C 2 + V 2 ) 2 

1 - \ 2 < 0. (38) 

4£ 2 (1 - k) 2 

Then we can draw a phase diagram of S-I transition with respect to rj and k (see Fig. 4). 



(|36D ~ (p8|) can be also expressed as follows 



1 / Ah 2 \ 

J xy = J z + 2 ( h + —j— ) @@ (critical region), @@@@ (39) 



1 / Ah 2 \ 

Jxy > Jz + - ( h -\ — — j @@ (superconducting region), (40) 



1 / Ah 2 \ 

Jxy < Jz + 2 [ h + @@( m sulating region). @@ (41) 
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The left-hand side is kinetic energy which causes quantum fluctuation. The right-hand side 
is potential energy which causes crystallization. In other words the former indicates wave 
property and the latter indicates particle property (localization) of Cooper pair, Therefore 
the phase diagram of our system is determined by the competition between them. 

5. SUMMARY 

We have focused on the S-I transition in ultra-thin films and considered randomness 
contribution as a factor to determine the border between the two phases (S and I). Su- 
perconducting order parameter (S x ) has been calculated as function of randomness rj, and 
influence of rj to S-I transition has been investigated. As a result, we have shown that (S x ) 
decreases with increasing i] and obtained analytically the critical value of randomness rj c . 

The existence of random potential causes scattering of Cooper pairs, and disturbs their 
free motion. Randomness rj becomes an important factor for the occurrence of electric 
resistance R. There is a one-to-one correspondence between the critical value of R and 
rj (R c and r) c ). Therefore solving i] c is equivalent to solving R c . It has been well known 
that the value of critical sheet resistance in ultrathin films is on the order of /i/4e 2 , but its 
exact coefficients has not been calculated yet. Our results may have some relevance to the 
determination of R c . 

We have obtained the results that randomness and two-body interaction causes localization 
and destroys superconductivity. These are capable of explaining, at least qualitatively, the 
main characteristic of S-I transition. In order to proceed to a more quantitative discussion, 
it is necessary to employ a more refined approximation in evaluating higher-order Green's 
functions. 



9 



ACKNOWLEDGMENTS 

We are very grateful to S. Saito, K. Sano, and B. H. Valtan for useful discussions and 
critical reading. 



10 



REFERENCES 

1. R. T. Scalettar, N. Trivedi, and C. Huscroft Phys. Rev. B59, 4364 (1999). 

2. N. Trivedi, R. T. Scalettar, and M. Randeria Phys. Rev. B54, R3756 (1996). 

3. M. P. A. Fisher, P. B. Weichman, G. Grinstein, and D. S. Fisher Phys. Rev. B40, 546 
(1989). 

4. P. Phillips, Nature, 406, 687 (2000). 

5. J. H. Sch6n, Ch. Kloc, and B. Batlogg, Nature, 406, 702 (2000). 

6. N. Markovic, C. Christiansen, and A. M. Goldman, Phys. Rev. Lett. 81, 5217 (1998). 

7. N. Markovic, C. Christiansen, A. M. Mack, W.H. Huber, and A. M. Goldman, Phys. 
Rev. B60, 4320 (1999). 

8. J. S. Penttila, U. Parts, P. J. Hakonen, A. J. Jordan, and R. P. Barber, Phys. Rev. 
B82, 1004 (1999). 

9. T. Yamaguchi, R. Yagi, A. Kanda, Y. Ootuka, and S. Kobayashi, Phys. Rev. lett. 85, 
1974 (2000). 

10. R. Fazio and G. Sch6n, Phys. Rev. B43, 5307 (1991). 

11. G. M. Zassenhaus and J. D. Reppy, Phys. Rev. Lett. 83, 4800 (1999). 

12. I. F. Herbut, Phys. Rev. B61, 14723 (2000). 

13. I. F. Herbut, Phys. Rev. B57, 13729 (1998). 

14. L. Zhang and M. Ma, Phys. Rev. B45, 4855 (1992). 

15. K. G. Singh and D. S. Rokhsar, Phys. Rev. B46, 3002 (1992). 

16. J. K. Freericks and H. Monien, Phys. Rev. B53, 2691 (1996). 

17. T. Matsubara and H. Matsuda, Prog, Theor. Phys. 16, 569 (1956). 

11 



18. H. Matsuda and T. Tsuneto, Suppl. Prog. Theor. Phys. 46, 411 (1970). 

19. K. S. Liu and M. E. Fisher, J. Low. Temp. Phys. 10, 655 (1973). 

20. D. N. Zubarev, Usp. Fiz. Nauk 71, 71 (1960). 

21. J. Kondo and K. Yamaji, Prog. Theor. Phys. 47, 807 (1972). 

22. H. B. Callen, Phys. Rev. 130, 890 (1963). 

23. PW. Anderson, Phys. Rev. 112, 1900 (1958). 



12 



FIGURES 

Fig. 1. Two-dimensional square lattice. The four sites i+a n (n = 1 ~ 4) are the nearest-neighbor 
sites of the site i on the square lattice. 

Fig. 2. The behavior of superconductivity (S x ) with increasing randomness rj. (S x ) continuously 
decreases and completely disappears at each specific point. Two body interaction k contributes to 
advancing the transition from superconducting to insulating. 

Fig. 3. The behavior of superconductivity (S x ) with increasing two body interaction k. (S x ) 
continuously decreased and completely disappeared at each specific point. Randomness rj con- 
tributes to advancing the transition from superconducting to insulating. 

Fig. 4. The phase diagram for our lattice boson model with randomness rj and two-body inter- 
action k. The borderline between two phases (S and I) is determined by the competition between 
wave property and particle property. When wave property is stronger than particle property, the 
system is in superconducting. When particle property is stronger than wave property, the system 
is in insulating. 
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